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Abstract 

We consider the problem of approximating two, possibly unrelated probability distributions from a single 
. complex-valued function and its Fourier transform. We show that this problem always has a solution within a 

specified degree of accuracy, provided the distributions satisfy the necessary regularity conditions. We describe 
the algorithm and construction of tp and provide examples of approximating several pairs of distributions using 
the algorithm. 

^ 1 Introduction 

In this paper we consider the problem of reconstructing two different probability distributions from a single complex- 



valued function with a specified degree of accuracy. This problem was suggested by the interpretation of the wave 
function ^ which is a solution of the Schrodinger equation. As is well known (x)'^ (x) is the probability density 
in a position space and ^>*{p)'^(j)) is the corresponding probability density in a momentum space, where ^'*(a;) is 
the complex conjugate of '^{x) and ^'(p) is the Fourier transform of Ij. We raise the related question as to 

^ whether any two distributions can be approximated within a specified error from a single function by calculating its 
' modulus or that of its Fourier transform, calculations like those used in Quantum Mechanics. 
I Our answer to this question is affirmative. 

' We prove this assertion by constructing a complex-valued function ^ which approximates two given distributions 

I within a given error, provided the two distributions satisfy some regularity requirements specified below. In particular, 
we show that the cumulative distributions can be approximated pointwise using the indefinite integral of the modulus 
of ■0 or of -0. 

The paper is organized as follows. In the next section we list assumptions, conventions and notation. In section 

3 we first introduce an expression for evaluating the error in the approximation of the two given distributions from 
the modulus of -0 or its Fourier transform. We then show how to construct this function. In the next section we list 

4 pairs of distributions approximated using construction described in section 3. In section 5 wc discuss the method, 
limitations and possible generalizations. In the appendix we sketch a proof that the function approximates any 
two given distributions on subintervals as claimed and a proof that the cumulative distributions can be approximated 
pointwise (uniformly on closed intervals) using the indefinite integral of the modulus of or of "0. 

2 Definitions, conventions and assumptions 
2.1 Assumptions 

1. X and P are any two variables; fx{x) is the probability density function for X with corresponding cumulative 
distribution Fx {x) ; fp (p) is the probability density function for P with corresponding cumulative distribution 
Fpip) 

2- fx{x) and fp{p) have finite variances and they (and their square roots) are continuously differentiable with 
Fourier transforms in and in (for definitions of and see 2 ; throughout, we denote the norm in 
aslllb) 

3. for every e > 0, there exist Xi, and Pt such that: i) for all x G [—Xb, Xb] and for all p e [—Pb, Pb] , fx{x) > 
and fp{p) > 0; and, ii) Fxi-Xb) < e/16, Fp(-Pb) < e/16, FxiXb) > 1 - e/16, and Fp(Pb) > 1 - e/16. 
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2.2 Definitions and Notations 

1. Let n be any integer, set J = 2", Ax — Xi,/J, and Ap = Pb/J, where Xf, and Pi, are given in assumption 3; if 
necessary, increase n (and/or Pf,) so that Fp{Ph — Ap/2) > 1 — e/16 ; 2J is the number of intervals of width 
Ax and Ap in [-Xf,, X^] and [-Pt - Ap/2, Pj, - Ap/2] 

2. s is a number defined in the Construction below; it indicates the number of subintervals into which each interval 
Aa; is divided. It is used to control the size of the error in the approximations 

3. t = s/{AxAp) 

4. j and k are integers, between — sJ, and sJ — 1; /, m and m' are integers between — J, and J — 1 or J, unless 
noted otherwise 

5. 4>s{x) is a function of the form used in the "sine approximation" [3] , [4]: 

sj-l 

{^) = X! V fx{jAx/s) sine 

j=-sj 

where 

, , sin(7ra;) 
smc(a;) = 



fx — jAx/s\ 
[ Ax/s 



6. Qx = max fx{x), for x e [-Xt, Xt]; Px = min fx{x) > 0, for a; G [-^b, -'^b]; ^'p = min fp{p) >0, for 
p e [-Pb - Ap/2, Pb]; Q'x = max dfx(.x)/dx, for a: G [-X,,, Xb] 

7- aj ^ f x{j Ax/ s)Y'^ Jot j^-sJ^-sJ +!,■■■ , s J - 1 

8- P^.;j,;. = a^, for < z < s • J - J 

9. prp;, = //ij'^^^ /p(z)dz for - J < ; < J - 1 

10. /(.) denotes the Fomicr transform of a function /(.), /*(.) denotes the complex conjugate of a function /(.) 

11. Mj is a function, defined recursively in the Construction, mapping the integers J G { — sJ, —sJ + 1, sJ — 1 
} to the numbers {—sJ/Ax, —s{J — 1)/Aa::, s{J — \)/ Ax, sJ/Ax} 

12. ruj — Mj ■ Ax/s; rrij is an integer 

13. Si = { i : i ^ {-sJ, -sJ + 1, sJ - 1} , and Mj = I ■ s/ Ax } for i = - J, ~J + 1, J. Equivalently, 5; = 
f^]^{l ■ s/ Ax), where f^l{ ) is the inverse mapping 



3 Construction 



Given fx{x), fp{p), £ and n, and X^ and Pb with the properties described in Assumptions we construct a complex- 
valued function i}} which satisfies the following two inequalities 



,7-1 

E 



t-^{qt)-r{qt) - fp{q)]dq 



< e 



(1) 



and 



j-i 

E 

i=-j 



(i+l)Aa: 



• - /x(a;) Ida; 



< £ 



(2) 



where i is a scale parameter defined as 
s 



t = 



AxAp 



(3) 
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and s is defined to be the smallest integer so that inequalities (H])-® are satisfied: 



Ax £ 



S 



1/2 



3 + 2-log{2sJ) 



< 



4(2J+1) 



||0,(x) - Vfxix)h < £/8(2J+ 1)4 



^-^JmJ^{2 + 2-log{2sJ)] < 



A{2j + iy 



(4) 

(5) 
(6) 

(7) 



fx{x)dx~ J2 fxijAx/s) — 



-Xb 
sJ-1 

E 



j = -sj 

(i+l)Ax/s A 

fx{x)dx^ J2 fx{jAx/s) — 



j = -sj 



< 



(8) 



Many of the functions already defined, such as (j)s{x), or to be defined, such as ip{x) and Mj, depend on s and n. 
For brevity we suppress explicit notation of that dependency. 

The next step in the construction of ip is to define subintervals. First we split the interval [^Xi,, Xh] into 2J 
intervals of length Ace: [I Ax, {I + l)Aa;], / = — J, — J +1, • • • , J — 1. Next we further divide each interval [/Ax, {I + 
l)Ax] into s subintervals and number them from j — — s J, — s J+ 1, to s J— 1. Define (j = — s J, — s J+1, • • • sJ— 1) 



Aa; 



X 



jAx 



1/2 



(9) 



We continue the construction by introducing the function Mj which maps the integers j € { —sJ, —sJ+ 1, s J— 1 } 
to the numbers {— sJ/Ax, — s( J— 1)/Ax, s(J— 1)/Ax, sJ/ Ax}, using the following recursive algorithm: initialize 
j to — sJ and / to — J. Then proceed with the following steps: 

I Set pTx-j.z = for z — — 1 and set prx-j,z — Ylim=j each integer 

z £ { 0, I ■ ■ ■ sJ — j }, and set 

Prp;i = lil-I^Ap fp(^)dz for - J < / < J - 1 

II If an integer z exists between —1 and sJ — j such that the inequalities < prp-i — pr^.j ;, < e/(2(2 J + 1)), and 

P'f'x-.j.z+i ~ Pfpd > hold, then set k ~ z; otherwise, set k — sJ ~ j 

III If fc > then set all Alj, Mj+k equal to Is/ Ax 

IV Set j ~ j + k + 1 and Z = Z + 1; if j = sJ, stop, otherwise go to step V 

V If Z + 1 = J, set Mj, Mj + 1, Msj^i all to sJ / Ax and stop; otherwise return to step I. 

Step II always returns a value of fc. The iterative algorithm continues until stopping, either at step IV or V. 

Properties of Mj that will be used subsequently include: for each j G —sJ,...,sJ— 1, the mapping Mj is a 
non-decreasing function of j; each Mj ■ s/ Ax is an integer (i.e. Mj ■ Ax/s = nij for some integer rrij € — J, J ). 

Using this definition of Adj, we can now define ip{x) appearing in equations (IT])-([2]) as: 



ip{x) = 

s.J-1 

E 




S.J-1 

E "j^" 

3 = -s.J 



TX 



jAx 



x — j Ax/s 
Ax/s 

X — j Ax/s 
Ax/s 



(10) 
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Apart from exp{2iTTMjx) and the sJth term in the series (ITUl) . the function iIj{x) represents the Whittaker-Shannon 
interpolation formula '5' -'6^ for •\/ fx (x) ■ 

In addition to approximations given in ([l])-© we also prove uniform approximations to the cumulative distribu- 
tions: let us define 

rP 

fp{q)dq 

-Pb-Ap/2 



Fp{p) 
Fp{p) 



Pb-Ap/2 



t ■ ilj{qt) ■ i}*{qt)dq 



(11) 



and 



Fx{x) 
Fx{x) 



fx{z)dz 
ijjiz) ■ ip*{z)dz 



(12) 



Then, it can be shown that with fx{x), fp{p), and e* given and with Xt and Pt having the properties described 
in Assumptions, one can choose n large enough that the complex- valued function given by the construction 
can be used to uniformly approximate the cumulative distributions Fx(x) and Fp(jp) for x G [— Xb,Xf,], and for 
[-Pb- Ap/2,P;,- Ap/2]: 



and 



Fx{x)-Fx{x) 



Fp{p)-Fp{p) 



< e* 



<e* 



(13) 
(14) 



For derivation of ([T]) and (|13p see Appendix. Approximations ([2]) and ()14p are derived in a similar fashion. 



4 Examples 



To illustrate the approximation, we compare the integral of the probability density function over a series of intervals 
with the correpsonding approximation generated by ip from ([Ij, ([2]) and ()10|) . Specifically, in P space we plot 

Ip-Ap/2 fpil)d<l and the corresponding approximation ^p'^^^j'^^ ■ '4'iQt)ip* {qt)dq. Similarly, in X space we plot 

Ix-Ax/2 fx{q)dq and the corresponding approximation ^^^'^^j'^ ' ^*il)dq] we approximate the last integral by 
ip{x)ili* {x)/Sx, as Ax is small. 

We applied the method to four pairs of distributions; for each pair, we use a single ip (and its Fourier transform) 
to approximate them with n=7, s—QA and Aa; = Ap — 5/2" k, 0.039: 

1. Figure 1 illustrates the approximation using Equation (jlOp. for the distributions 



fx{x) 
and fpijp) 



/2ti 



0.75e-«(^+5f +o.25e-8(-i)^ 
1 



'2tt 



2. Figure 2 illustrates the approximation using Equation (|10p . for the distributions 



fx{x) = e ^, a; > and Jp{p) = 



1 



3. Figure 3 illustrates the approximation using Equation (|10p . for the distributions 
1 



fx{x) 



/27r 



and /p(p) = e-P, p > 



4. Figure 4 illustrates the approximation using Equation (|10p . for the distributions 

/2 



fxix) = and fp{p) = -^e'^' 
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5 Discussion 



Based on the algorithm described in section 3 we have shown how to approximate two given distributions from a 
single complex- valued function with a specified accuracy. The main result of our paper is given in expressions , 
P3m4p and pU)) . Expression ^TU\\ can be viewed as an extension of well-known numerical methods based on the sine 
function providing a method for approximating two distributions simultaneously. 

The construction and arguments given here show that there always exists a function of the form given by pop 
which approximates two given distributions with a specified degree of accuracy provided our assumptions (1-3) are 
satisfied. For successively higher degrees of accuracy, the method allows construction of a sequence of approximating 
functions but this sequence may not necessarily converge, say in L^, to a limiting function. 

A modification of our approach or a similar method based not on the sine function but on some other basis 
or wavelet series [7] might be used to approximate the two distributions and converge to a limiting function. For 
example, one might attempt to use the Hermite polynomials multiplied by the appropriate Gaussian function as a 
basis. 

It remains speculative, however, whether our approach can be modified so as to approximate arbitrary pairs of 
distributions with a basis other than the sine functions. Further speculation about this potential generalization is 
beyond the scope of this paper. We would like to stress that the main goal of the present work is not to find the 
best numeric approximation (though we can always meet an increasing demand in accuracy) but to establish the 
existence of an algorithm allowing uniform approximation of the cumulative distributions. 



APPENDIX 

We sketch the derivation of ^ and (|13p . First we write down the Fourier transform of ip{x): 
^(g) = J V(x)e-*-2'^-^-«(ia; 

sj-l 
j=-sj 

Ax/ s 




Ax , , , x\ 



e 



-i-27r-jA£(9-Af,) 



(A-1) 



In (|A-ip i? is a rectangular function: 



{1 when — i < g < i, 

i wheng = ±i, (A-2) 
otherwise 

We consider 

= / t-^{qt)-riqt)dq, (A-3) 

the integral of t ■ ipiqt) ■ 'ip*{qt) over the interval [{I — ^)Ap, (/ + ^)^p], for each I e — J, — J + 1, J — 1. Using 

expression (jA-ip . rescaling the integration variable, and using definition of integer rrij : rrij — AxMj/s, for this 
integral we obtain: 



(,+ 1) sJ-l 

I{1) = dq' } ajakR{q' - m,)R{q' - mk) 



r('+5) ^izi 
where aj — [Axfx{i^x/s)/sY^'^. 
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Property (IA-2[) of the rectangular function R implies that each j, k term in the integrand is either for the full 
range of integration or vanishes for all but a subinterval of length 1. The latter holds only if rrij = mk and we have: 
I < rrij ± i < {I + 1) if and only if / ~ nij, for some integer rrij. After interchanging the order of integration and 
summation, each integral in the sum is due to the factor e^''^^'^^^^''''^ , unless j ~ k. On the other hand, if j = fc 
the integral is just a|. Thus, we can write (I — — J, — J + 1, J — 1): 



/(/) 



(i + i)Ap 



(A-4) 



where we recall the definition of Si: Si is the set of all j S —sJ,~'sJ+ 1, sJ for which Mj ■ (Ax/s) — I. 

There are 2 J terms /(/) of the form given in (E^ : the approximation given by inequality ([T]) will be proven if 
we can show that for j = — J, — J + 1, J — 1 the following relation holds: 



jes, 



E 

3<^Si 



a? - / fp{z)dz 



< 



2J 



(A-5) 



Either inequality (jA-Sp holds for I = —J,—J+ 1, J — 1, and we are done, or there is a smallest I violating (jA-Sp . 
say h, such that Ejes,, l^j - I{i^^}Ap /p(^)^^l 

> e/{2J). This implies that ji > sJ — 1, where ji = max{j : j € Si]. Indeed, if this were not the case (i.e. 
ji < sJ—1), than Afjj+i would have been defined by the algorithm described in section 3 as li since a^j+i < e/ (2J+1) 
and ji + 1 would also have been in Si-^ . This would be a contradiction since ji was supposed to be the largest element 
of Si^. Therefore, ji > sJ — 1. 

By choice of s in section 3, requirement dH): | X^i^-j Sjes, '^j ^ IjAx^^f'^^'^ fx{x)dx\ is less or equal e/8, and 
because fx{x) integrates over the interval [— Xfc,Xf,] to at least 1 — e/16 (Assumption 3), we have: 



^ a]>l- e/16 - e/8 > 1 - e/4 



(A-6) 



Since fp{p) is a probability density function and by Assumption 3 and the algorithm described in section 3: 



1 > / fp{z)dz = V / fp{z)dz 

J-P,-lAp i^-jHi+h^P 

i=~JjeSi j=-sj 



(A-7) 



The last inequality above is simply inequality (jA-6p . By the Construction, jl^l~^i^^p fp{z)dz > X^jgs, ^'^^ every I. 



This fact and (|A-7p imply: 



fpiz)dz - E E = 



E 

i=-j 



fp{z)dz 



jes, 



1 



(A-8) 



Further, the upper bound of 1 in expression (|A-6p implies that: 
Jp{z)dz = E / fp{z)dz 

(ii + i)Ap i^i^^^J{l-^)Ap 

3 



(A-9) 
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Now Si is empty for / > li since Mgj-i = h and Mj is non-decreasing. In other words X^jes, a| = for / > li. 
This fact and the inequahty in Expression (jA-91) yield: 



,7-1 



J-1 .(/+i)Ap 

V / fpiz)dz-0 



3 



(A-10) 



Finally, inequalities (|A-8I) and (lA-lOp lead to the desired result, proving inequality ([T|). 

Proof of inequality ([2]) requires no new qualitative features but straightforward tedious calculations; we do not 
cite it here. 

To prove the uniform approximation for the cumulative distributions, inequality (|13p. choose n large enough 
that Qx ■ Ax = Qx ■ A^b/S" < £*/2, construct '4){x) so that inequalities ((T)) and (0) hold with £—£*/2 and let 
5{x) = Fx{x) ~ Fx{x) . We now show that 5{x) < e* for all x £ [jAx, (j + 1)A.t] and for j = - J, -J + 1, • • • , J- 



1. Suppose the maximum of 



Fx{x) — Fx{x) occurs at X = xq, for xo G [jAx, (j + l)Aa;]. Let Mp^ denote 



Max{Fx{x),x e [jAx, {j + l)Ax]}. If Fx(a;o) > -Fx(a;o), then 
S{x) <Mf^ -FxijAx) 



(A-11) 



since Fx{x) is non-decreasing. Further, the mean value theorem implies Mp. < Fx (jAx) + Qx ■ Ax since Qx is the 
maximum of fx{x). Using this result in (jA-ll[) gives 



S{x) < Fx{jAx)+Qx-Ax~Fx{jAx) 

< e*/2 + e*/2^e* 

since Fx{jAx) - Fx{jAx) < e*/2 follows from ^ and IQ^^ • Aa;| < e*/2 by choice of 
On the other hand, if Fx{xa) < Fx{xo), then 

S{x) < Fx{{j + l)Ax)~Fx{3Ax) 

< FxiU + l)Ax)+e*/2-FxijAx) 

< Fx{jAx)+Qx- Ax + e~ Fx{]Ax) 

< e*/2 + e*/2^e* 



(A-12) 



(A-13) 



The chain of estimates above follows from the fact that Fx{x) is non-decreasing, Qx is the maximum of fx{x), and 
by choice of n. 



The proof for 



Fpip)-Fpip) 



is similar. 
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Figure 1: Approximation to the integrals /jfj!^jf/2^ fp{l)d(l (top) and J^^^^j2 fx{q)dq (bottom) 
based on Equation pOjl for the pair of distributions fp{p) — exp(— p^)/\/27r, and fx{x) = 




-2-1012 -2-1012 

P X 



Figure 2: Approximation to the integrals fpil)d<l (top) and ^^^^^^2 fx{<l)dq (bottom) 

based on Equation pup for the pair of distributions fp{p) = '^4^6"'' ^ fx{x) = for 
a; > and fx{^) = for a; < 0. 



Figure 3: Approximation to the integrals Jpl!^^pl2 fpil)'^1 (bottom) and ^^^^^^2 fx{q)dq (top) 
based on Equation pO|) for the pair of distributions fx{x) — --^e~^ and fp{p) ~ exp(—p) 
for p > and fp(p) ~ for p < 0. 



Figure 4: Approximation to the integrals J^l^^^pj^j fp{l)'^1 (bottom) and j^^^^j2 fx(l)'^1 (^op) 
based on Equation (|10p for the pair of distributions fp(v) = —^e~^ ^'',a,nAfx(x) — —^e~^ . 



